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Abstract. A topological constraint on the dynamics of a magnetic field in a 
O j' flux tube arises from the fixed point indices of its field line mapping. This can 

explain unexpected behaviour in recent resistive-magnetohydrodynamic simulations 
of magnetic relaxation. Here we present the theory for a general periodic flux tube, 
■ representing, for example, a toroidal confinement device or a solar coronal loop. We 



show how an ideal dynamics on the side boundary of the tube implies that the sum 
of indices over all interior fixed points is invariant. This constraint applies to any 
O ■ continuous evolution inside the tube, which may be turbulent and/or dissipative. We 

c/3 \ also consider the analogous invariants obtained from periodic points (fixed points of 

the iterated mapping). Although there is a countably infinite family of invariants, we 
show that they lead to at most two independent dynamical constraints. The second 
constraint applies only in certain magnetic configurations. Several examples illustrate 
the theory. 
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1. Introduction 



X 

Magnetic flux tubes are basic structural features in astrophysical magnetic fields like 
the Sun's atmosphere, or may represent the toroidal magnetic field in thermonuclear 
confinement devices such as tokamaks pQ. Within a single flux tube, the magnetic 
sub-structure may be extremely complex, particularly when the plasma undergoes a 
turbulent evolution. A useful technique for characterizing the magnetic field in such a 
tube is the field line mapping from one end to the other (defined in section [2]). This 
reduces the field to a two-dimensional function that is simpler to analyze, yet retains 
much important information. Field line mappings have been used in solar physics, 
where strong mapping gradients are a likely indicator for preferred locations of magnetic 
reconnection [2], and to model laboratory devices with toroidal geometry (3J H]. The 
latter have exploited the mathematical interpretation of a flux tube as a dynamical 
system in two space dimensions, with the role of "time" played by the direction along 
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Figure 1. Initial and relaxed magnetic fields in the "braided magnetic field" simulation 
(see [Till E] for details and background information) . Only selected magnetic field lines 
are shown (the field fills the whole box), traced from the same starting points on the 
lower boundary in each case. 

the tube. Powerful techniques from dynamical systems theory may thus be applied to 
study the three-dimensional magnetic structure at a given time snapshot. In this paper, 
we show how a basic topological property of the field line mapping — the fixed point 
index — may lead to a physical constraint on the time evolution of a magnetic flux tube. 

This topological constraint applies under any evolution where field lines on the 
side boundary of the flux tube (which is by definition a magnetic surface) undergo an 
ideal evolution that preserves their connectivity [5j. Essentially, the interior dynamics 
are restricted because the total index of interior fixed points is determined by that 
of fixed points on the boundary (section [3]). The constraint does not depend on 
the nature of the dynamics in the tube, which may follow the equations of ideal- 
magnetohydrodynamics (MHD), resistive MHD, or any other continuous evolution of the 
magnetic field. However, the work was initially motivated by resistive-MHD simulations 
of magnetic loops in the solar corona (atmosphere). The coronal magnetic field builds 
up stress and energy as it is twisted by motion of its footpoints in the solar interior. 
The key question is: what will be the final state of its subsequent turbulent relaxation? 
Knowing the appropriate constraints that apply during this relaxation would allow us 
to predict the final state, and hence place limits on the amount of magnetic energy that 
can be dissipated during the relaxation. This has direct implications for understanding 
how the magnetic field heats the corona to extreme temperatures [7j. 

The appropriate constraints for a resistive-MHD relaxation remain under debate [8] . 
The theory of Taylor [9], which predicts relaxation to a so-called linear force- free field, 
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Figure 2. The physical setup: a magnetic field B on a cylinder of radius R. 



can explain the final state in certain laboratory experiments, and has been conjectured 
to apply to the solar corona [TJ. However, its applicability to astrophysical magnetic 
fields like the corona is unclear, with recent numerical simulations of a solar coronal loop 
[TO} [TT] finding a final state in conflict with expectations from Taylor theory. Figure [I] 
shows one of these simulations, which we refer to as the "braided magnetic field" (due 
to the initial structure) and use as an illustrative example throughout this paper. We 
have recently shown that the index constraint can explain why these simulations are 
unable to reach the predicted Taylor state [T2] . 

In this paper, we develop fully the initial idea of [12] for a general flux tube, 
explicitly accounting for the (important) influence of the side boundary and exploring 
the role of higher periodic points. The latter are fixed points of iterations of the field 
line mapping. To allow for these, we focus on a flux tube with periodic boundary 
conditions (i.e., the distribution of the normal magnetic field is the same across both 
ends of the flux tube). In that case, the tube end is a Poincare section in the dynamical 
system analogy, and the field line mapping is a return map. Periodicity is natural if 
the cylinder represents a toroidal geometry, but is also typical in numerical simulations 
of solar coronal magnetic loops where the initial conditions are given by an analytical 
equilibrium and the ends remain "line-tied" throughout the evolution. In fact, we find 
that in most cases the higher periodic points do not impose any additional constraint 
on the dynamics. This result is formulated in Theorem 14.21 and proved in section HI 

2. The magnetic field line mapping 

We shall consider magnetic fields B(x) = B r e r + B^e^ + B z e z on the cylinder 
{0 < r < R, 0, < < 2tc,0 < z < 1}, sketched in figure [2J The lower and upper 
boundary discs are denoted D and D\ respectively. We impose the following conditions: 



B^O, 

B r \ r=R = 0, 

B z \d — B z \d x > 0. 
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In sections E] and H] we will consider time-dependent magnetic fields B(x, t), and these 
conditions will be applied at all times. In particular, ffTcj) allows the magnetic field on 
D or Di to be time- dependent, providing that B z remains periodic. Static "line-tying" 
is a special case. Since V • B = 0, the magnetic flux is the same through any horizontal 
cross-section of the cylinder, which therefore represents a magnetic flux tube. Moreover, 
all magnetic field lines (integral curves of B) connect D and D\ in the same direction, 
and can be parametrized using the coordinate z. The field line starting at x G D is 
denoted F(x , z) and found by integrating 

dF(x ,z) = B(F(xp,z)) 

dz S 2 (F(x ,z)) 1 ] 

from z = to z = 1, with F(xo,0) = xo. Tracing field lines from every xo yields a 
mapping / : D — > D x given by 

/(xq) :=F(x ,1), (3) 

called the field line mapping of B. Since we are interested in physical magnetic fields 
with finite dissipation, we assume that B is differentiable, and hence that / is a 
diffeomorphism of the disc to itself. It follows from the fact that / is a homeomorphism 
that the boundary circle 8Dq is mapped to itself. 

From V • B = 0, it may be shown that the Jacobian determinant of / is 

where Df^ := df 1 jdxK Condition ( TTcl) then implies that \Df \ > 0, so the mapping 
is orientation-preserving. In those fields where _B z (/(xo)) = 5 z (xo) for every field line, 
/ is also area-preserving. Such symplectic maps have been widely studied as models 
for magnetic fields in the fusion context [13J, but a general field line mapping is only 
symplectic when written in appropriate canonical coordinates [14] . Although most of 
our examples do retain this property in physical coordinates, the results in this paper 
do not depend on it, requiring only conditions (iTaj) - (iTcj) . 



3. Fixed points 

To analyse the local behaviour of the field line mapping / we introduce, on the disc D 1} 
the mapping 

V/(xq) = /(x ) -Xq. (5) 

Although V/ is not a proper vector field, we can treat it like one for the purposes of this 
paper. The main use of vj will be to define the fixed point index, but we note first that 
it leads to a simple but effective "colour map" technique for visualising the 2-d mapping 
/ [15]. In Cartesian coordinates, each point x G D is assigned one of four colours, 
according to the two Cartesian components of v^: 

red if v x f > and v y f > 0, 
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Figure 3. The uniform twist field (j6]) with (a) a = 7r/4 and (b) a = 7r. The vectors 
show the direction of V/, normalised to unit length for clarity. Here R = 2. 



yellow if Vf < and > 0, 
green if < and ujl < 0, 
blue if > and ^ < 0. 

As a simple example, figure [3] shows the "uniform twist" magnetic field 

B = B twist + e z , B twist = -aye x + axe y , (6) 

for two values of the twist-angle parameter a. The arrows in figure [3] show the direction 
of Vf, which for this field may be shown to be 

Vf = (cosa)xo — (sma)yo — xo, (7a) 
v V f = (sina:):Eo + (cosa)?/o — Vo- (7b) 

A more complex example is given by the braided magnetic field simulation 
introduced in section [TJ A sequence of colour maps obtained at different times during 
this simulation is shown in figure HI It neatly reveals how the field line mapping is 
simplified during the relaxation. 

Central to our analysis are the fixed points where /(x ) = x , or equivalently 
Vf = 0. Isolated fixed points show in the colour map as intersections of the curves 
= and v y f = 0, i.e., points where all four colours meet. Indeed, the colour map was 
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Figure 4. Time sequence of colour maps from the braided magnetic field simulation 
(see [IT] for details of the simulation). Letters y, g, b, r indicate the colours yellow, 
green, blue, and red (for readers viewing in greyscale). 



introduced to aid the numerical search for fixed points pj)]. The uniform twist field (J6j) 
has a single fixed point at the origin x = y = 0, while the braided magnetic field has 
multiple fixed points at different times. 

The reason for focusing on fixed points of / is simple: they are topological properties 
of the mapping, and hence of B. In an ideal time evolution, the fixed points are 
invariants. But even under a non-ideal evolution they cannot be arbitrarily created 
or destroyed. Importantly, we can make a stronger statement not just about the fixed 
points in isolation but about the local structure of the mapping around them. This 
is expressed through a topological property called the (Poincare/-Hopf) index of each 
fixed point. 



3.1. Index of a fixed point 

The index of a fixed point of / depends on the local mapping around the fixed point, 
so to assign an index to fixed points lying on the boundary dD , the mapping / must 
be extended outside the disc r = R. To this end, following [16], define the extended 
mapping 

/(*,) = ( /(X0) ifW£fl ' (8) 

I /(iWW) if |xo| > R, 
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which is continuous and introduces no additional fixed points. (An alternative extension 
is used in [T7].) 

Let T be a closed curve enclosing an isolated fixed point x G D but enclosing no 
other fixed point. The index of x , denoted ind XQ /, is the winding number of the closed 
curve Vjp(r) about the origin [16]. It is an integer and so invariant under homotopy: it 
is a topological property of the local mapping around the fixed point. The extension 
/ means that this definition applies both to interior fixed points and to those on dD . 
Note that there are a number of equivalent definitions of the index |18j . The winding 
number may be expressed in integral form as 



ind 



2ir 



arctan | — 




1 

2^ 



vj dv y ~ 



v y - dv°~ 



(9) 



/ 



which is the two-dimensional case of the Kronecker integral [15]. Evaluated around a 
closed curve enclosing more than one fixed point, the integral yields the algebraic sum 
of their indices. 

The fixed point in the uniform twist example has index +1. As another simple 
example, figure |5] shows the "hyperbolic" magnetic field 

B = B hyp (r, 0) + e z , (10a) 



B hy P( 



2 cos(20) sin 



/7rr\ 

sm It) 



ixr 



■ cos 



sin(20)e^ 



(106) 



"KT\~ 
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comprising a uniform vertical component B z = 1 and horizontal components given by 
B hyp . The extension of the colour map outside the disc R = 2 (using /) is also shown. 
Here / has five interior fixed points and four on the boundary. The fixed point at 
^o = 2/o = has index —1 (hence the name "hyperbolic"), while the remaining four 
interior fixed points all have index +1. Of the boundary fixed points, the pair at xq = 
have index 0, while those at yo = have index —1. 



3.2. Global constraint 

The utility of the fixed point index is encapsulated in the following global result 
connecting the indices of all of the fixed points of / [19j [16]. We denote the set of 
fixed points of / by Fix(/). 

Theorem 3.1 (Hopf). If M is the closure of an open, connected subset of the plane 
whose boundary is a union of smooth simple closed curves and f, g : M — > M are 
homotopic maps with both Fix(f) and Fix(g) finite, then 

^2 ind vf = ind p9- 

In our case, M is the disc Dq, and the following simple argument shows that the 
value of the index sum for any / is unity [16]. Since the disc is contractible, all self-maps 
are homotopic. Thus any map / is homotopic to the "constant" map c : D — > D where 
c(x ) = for all x G D. This map has a single fixed point at x = with index 1. 
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Figure 5. Colour map of the "hyperbolic" magnetic field (|10a|) . with projected field 
lines showing the horizontal component B hyp (not the vector field v/). Black circles 
identify the nine fixed points of the field line mapping. 



Hence, by Theorem 13. 1[ the map / must also have index sum 1, provided that it has a 
finite set of fixed points. For our further analysis it will be useful to separate the sum 
of interior fixed points from those on the boundary. We define 

T^t = md p/' T d= ^2 md p/' ( n ) 

P eFix(/)nint(A)) peFix(/)naD 

giving us the following result. 

Corollary 3.2. If f is a continuous self-map of the disc D with Fix(f) finite, then 

T mt + T d = l. (12) 

In the absence of boundary fixed points, the interior sum T int is known as the 
Lefschetz number, or topological degree of the mapping /. For the hyperbolic field of 
figure |5l it is readily determined that T int = 3 and Tq = —2, verifying Corollary 13.21 for 
this particular example. 

The relation (fl2|) will be satisfied by the field line mapping of any continuous 
magnetic field B, at all times during its dynamical evolution. In general, a magnetic 
field will be able to exchange contributions between the two terms Tj nt and Tq, in such 
a way that their sum remains unity. However, if the physical situation is such that Tq 
remains fixed, then Corollary 13.21 implies that Tj nt must also be a conserved quantity 
in the dynamical evolution, irrespective of the details of the dynamics. This is not 
an unusual situation: for instance, in astrophysical plasmas — where the evolution is 
commonly ideal in most of the domain — it is often possible to choose a flux tube such 
that the boundary remains in the ideal region while the interesting (non-ideal) dynamics 
occur in the interior. 

The braided magnetic field simulation is a case where this topological constraint 
is significant and determines the final state of the system [12]. The sequence of colour 
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maps in figure S] shows how, although the number of fixed points decreases from 27 
at t — to 2 at t — 290 (through bifurcations which preserve the total index), the 
overall value T- mt = 2 is conserved. Note that, in this simulation, the evolution at the 
side boundary of the computational box remains ideal, so that the field line mapping is 
preserved there. Figure H] shows only part of the domain between x = ±3,y = ±3, but 
this region contains all of the interior fixed points and all of the non-ideal dynamics. 
The total index T int may be determined either by calculating the Kronecker integral (Q 
around the boundary of this square region [T5j, or simply by inspecting the sequence 
of colours around this boundary, which suffice to determine the winding number of v;. 
The latter is clearly seen to remain the same throughout the simulation. 

3. 3. Resistive diffusion: an apparent paradox 

For a fixed normal flux distribution on the boundary of the domain, it is well known 
that under resistive diffusion in a static medium, where 

f - < 13 > 

every initial field B must decay asymptotically to a unique minimum-energy potential 
field (satisfying V x B = 0). But, for a given normal flux distribution on the boundary, 
one can have fields with a wide variety of T int . How can T int be conserved if it differs 
from Ti n t of the minimum-energy field? 

The resolution of this paradox is suggested by the braiding simulation. Here the 
dynamical relaxation leads to a final state that preserves T; nt = 2, even though the 
minimum-energy field is B pot = e 2 , a uniform vertical field where the index is undefined 
(since all field lines are periodic). If the simulation is continued further, the horizontal 
magnetic field continues to decay, but only on the much slower global resistive timescale 
r = L 2 /r], where L is the global domain size. This decay suggests that the vertical field 
Bp t will be reached, but only asymptotically as t — > oo. This infinite limit is usually 
not of physical interest; the question is how the system evolves on the much shorter 
dynamical timescale. At all finite times, while any horizontal magnetic field remains, 
we expect T int to be preserved (providing that Tg is preserved). Further investigation of 
this point is merited in more general configurations. 

3.4- Generic fixed points 

What are the possible values of ind xo / for a fixed point? Non-degenerate interior fixed 
points correspond to non-degenerate critical points of V/, so have index either +1 or 
— 1. Non-degenerate means that det(Df — Id) ^ 0. A degenerate fixed point can always 
be decomposed by a small perturbation of the mapping into several non-degenerate 
fixed points, such that their total index equals that of the degenerate fixed point. Such 
degenerate fixed points are thus topologically unstable. In our physical system, we 
assume that all interior fixed points will be non-degenerate, or "generic" . 



Dynamical constraints from field line topology in magnetic flux tubes 10 




ind Xo / = ind Xo / = -l 



Figure 6. The local vector field for the two generic types of boundary fixed point, 
in the case that \Df\ — I at the fixed point. Dashed lines indicate the extension 
outside _Dq, which is determined completely by the mapping on 8Dq. 



The boundary condition ( |T5|) restricts the possible types of fixed point that we can 
have on the boundary dD . 

Lemma 3.3. Let x be an isolated, non-degenerate fixed point of f located on the 
boundary 8Dq. If \Df\ (xq) = 1 then ind yiQ f is either or —1. 



The proof is given in Appendix A , while the two possible cases are shown in figure EJ 
Observe that our periodic condition (flcj) . along with (TJJ, ensures \Df\ (x ) = 1 at every 
fixed point xo. Thus Lemma 13.31 applies to every fixed point on 8Dq, implying that 
Tg < for any /. So, by Corollary 13. 2 [ a periodic flux tube must have 

T in t > I- (14) 



4. Periodic points 

The success of Corollary 13.21 in explaining the final state of the braided magnetic field 
[T2] leads us to ask whether additional topological constraints arise from higher periodic 
points. This requires condition ( ITcj) that the flux tube is periodic. 

Writing f m for fofo...of iterated m times, a periodic point xo G -Do of (minimum) 
period m satisfies 

f m (x ) = x , (15a) 
/ n (x ) 7^ x for n < m. (15b) 

A periodic point is thus a fixed point of f m . In particular, the fixed points of / are 
periodic points with period 1. Figure [7] is a sketch of the second iteration, f 2 , for 
a hypothetical flux tube. Notice that the fixed points of f 2 include periodic points of 
periods 1 and 2. In general, the fixed points of a given iteration f q could include periodic 
points of any period m > 1 that divides q (written m\q). Notice also that the periodic 
points of period m fall into groups of size m corresponding to the same "periodic orbit" . 
In figured the blue periodic orbit comprises two period 2 points, xJj; and x^. 
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Figure 7. Field lines corresponding to some fixed points of / 2 . The point xi (in red) 
is a period 1 point (i.e. a fixed point of /). The points Xj, x^ (in blue) are period 2 
points, both belonging to the same periodic orbit. 



For a given iteration q we can consider (as before) the sum over all interior fixed 
points, 

P eFix(/9)nint(D ) 

and the analogous sum Tq over the boundary fixed points. 

Corollary 4.1. If f is a continuous self-map of the unit disc D , with Fix(f q ) finite 
for every q G N, then 

T int + T$ = 1 for all q G N. (17) 

If all Tq are conserved during the dynamical evolution, then the T q nt will be a 
countably infinite family of topological constraints on the magnetic field evolution. The 
question of interest physically is whether these additional constraints are important in 
practice, or, in other words, whether they are independent of one another. The following 
result shows that, in a large class of magnetic fields, they are not. 

Theorem 4.2. Let f be an orientation-preserving homeomorphism of the unit disc Dq, 
with Fix(f q ) finite for every q G N, and with ind Xo f q G { — 1, 0} for every periodic point 
xo on 8Dq. Then 

(i) If T} nt > 1 then T q nt = T} nt for all q G N. 

(ii) If T} nt = 1 then T-* nt = IMq G N if and only if the rotation number of f\dD 
irrational. 

According to Theorem 14.21 the index sums of higher iterations can differ from that 
of the initial mapping only if the initial sum is 1, and then only if the restriction of / 
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to the boundary dD has a particular form (namely a rational rotation number: this 
will be defined below). Thus for any other initial magnetic field, it will hold at all times 
throughout the ensuing dynamical evolution that T? nt = T^ nt for all q G N. This means 
that the higher constraints of Corollary 14.11 will be automatically maintained for any 
possible dynamics and thus pose no additional constraint on the evolution beyond that 
imposed by T^. 



4-1. Proof of Theorem 4-2 



The strategy is to use the fact that T? nt — 1 — Tg and focus on the boundary mapping 
f\dD - Firstly, we consider the sequence of indices of an individual fixed point under 
iteration of /. 

A fixed point x of / recurs as a fixed point of each iteration f q for q G N, and the 
index of this fixed point under each f q generates a sequence {ind X0 / 9 }. The possible 
index sequences for different classes of mapping have been quite well studied (see, e.g., 
[20J). If / is continuous and x is an isolated fixed point for each q, then Dold [2Tj 
proved that the sequence of indices satisfies the sequence of congruences 

Ufa) ind x f q/n = mod q (18) 

n\q 

for each q > 1 (called Dold relations). Here fi(n) is the Mobius function from number 
theory, defined by 

{0 if p 2 \n for some prime p, 

1 if n is square- free with an even number of prime factors, (19) 

—1 if n is square- free with an odd number of prime factors. 

Since the extended mapping / is continuous, we may use the Dold relations (|T8l) to 
obtain the following Lemma for fixed points on dD . Note that it follows from ([S]) that 
the fixed points of f q are the same as those of f q for all q. 

Lemma 4.3. Let f be a field line mapping whose boundary fixed points on each iteration 
have index either or —1. Then for each such fixed point x 0; 

indstof 9 = indxgf for all q G N. 



Proof. Induction on q (cf. a similar argument used by [22] )• Clearly the result holds for 
q = 1. Suppose that ind XQ / m = ind Xo / for all m < q. Now, separating the first term in 
the qth Dold's relation (II 8p gives 

/i(l) ind Xo r + ™<f q/n = mod q. 



n\q 
n>l 



Now /i(l) = 1, and since q/n < q the inductive assumption gives 
ind X0 / 9 + ind xo / ^ M™) = mod Q- 



n\q 
n>l 
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But 

n\q n\q 
n>l 

using a well-known property of the Mobius function. Thus we must have 

md xo f q = ind xo / mod q, 
and the result follows from the restriction that ind Xo / 9 G { — 1, 0}. □ 

Now consider the restriction f\dD of / to the boundary. This is an orientation- 
preserving homeomorphism of the circle. (Orientation-preserving here means that the 
order of any three points on dD cannot be permuted by /. Physically, it follows from 
the fact that / maps dD to itself, and the restriction that magnetic field lines cannot 
intersect when B ^ 0.) But an orientation-preserving circle homeomorphism can have 
periodic points of at most a single (minimal) period |23j. Thus there are two cases. 

Case (i): T^ nt > 1. In this case we must have Tq ^ 0, so that / has at least one fixed 
point on the boundary. But then / can have no boundary periodic points of any higher 
period: the fixed points of f q will be the same as those of / for all q G N. Applying 
Lemma 14.31 to each boundary fixed point gives T q nt = T^ t for all q. 

Case (ii): T? t = 1. This implies that / has no fixed points on ODq. Therefore 
there could be periodic points of some unique period m > 1. According to the theory 
of circle homeomorphisms, the existence of such points depends on the rotation number 
of the homeomorphism [23]. The rotation number of f\dD is defined as 

,. G k ((f> ) - 00 
P = llHl , 

fc-5>oo 2irk 

where <fto G dD$ and G is a lift of f\dD to the real line. The limit exists 
and is independent of the point (f) . Since f\oD is an orientation-preserving circle 
homeomorphism, p is rational if and only if f\dD has at least one periodic point. In 
that case, and if p written in its lowest terms is n/m, then all periodic points must have 
period m. Otherwise, p is irrational and f\oD nas no periodic points. We now apply 
Lemma [4.31 In the irrational case, we have Tq — for all q G N. In the rational case, 
we have 

n = { Tm ifm|g ' (20) 

I otherwise, 

where T m < is a constant given by the index sum of the mth iteration (the first with 
fixed points). This completes the proof of Theorem 14.21 □ 

4-2. Examples 

The following examples illustrate the two possible cases in Theorem 14.21 Case (i) 
where T^ nt > 1 is exemplified by the braided magnetic field simulation discussed earlier. 
Figure [8] shows computed colour maps for /, f 2 , and f 3 in both the initial and final 
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(a) initial state 




(b) Final state 




Figure 8. Colour maps for the first three iterations of / in (a) the initial state and 
(b) the final state of the original braided magnetic field simulation. 



states of this simulation. At both times, T^ nt = 2, and we see that T-^. — TJnt — 2, in 
accordance with Theorem 14.21 Thus the higher T? nt are automatically preserved in this 
simulation because T^ t is preserved and the dynamics are ideal near to 8Dq. Notice 
that, although the index sums are equal for each iteration, this does not preclude the 
appearance of more periodic points at each iteration. This is visible particularly in the 
final state, where there are 2 fixed (period 1) points, 16 period 2 points, and 36 period 
3 points. However, the additional periodic points appear in equal numbers with index 
+ 1 and —1. 

To illustrate case (ii) in Theorem 14. 2[ we begin by constructing a magnetic field 
whose field line mapping has T^ nt = 1 but T^ t ^ 1. This may be done by taking the 
uniform twist field B twlst + e z and adding a small perturbation. We choose a = it so 
that every point of D is a period 2 point of B twlst + e z . The perturbation takes the form 
of B hyp , except that it is rotated as z increases at the same rate as the uniform twist. 
This destroys the symmetry of the uniform twist map, leaving a finite set of period 2 
points, which correspond to the fixed points of the original B hyp field. Our overall field 
is 

B = B twist + eB hyp (r, - irz) + e z , (21) 

where e = 0.1 and a = n. The boundary mapping has rotation number p = 1/2, so has 
periodic points of minimal period 2 only. This is verified for the first four iterations of 
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Figure 9. Sequence of colour maps for iterations of the perturbed uniform twist 
magnetic field (j2"Tj) . where all four periodic points on the boundary have period 2. 



/ by the colour maps shown in figure |9j In particular, we find T? nt = 3, and so for this 
mapping 

| 1 for q odd, / 

TL={^ , (22) 
I 3 for q even. 

If a magnetic field with this boundary mapping were to evolve subject to ideal 
dynamics at the boundary, then our theory predicts two constraints on the evolution. 
We have verified this prediction through a new numerical simulation, solving the same 
resistive-MHD equations as the earlier braiding simulation (T0J ITT] . To initiate a 
dynamical relaxation, we add to the background field (12T|) six twist regions. These 
twist regions match those in the initial condition of the original braided magnetic field 
simulation (see [24J for details). Their effect is to add significant complexity to the field 
line mapping, triggering the formation of thin current sheets in the relaxation process. 
(Since the initial magnetic field is smooth, it must remain smooth for any finite time 
t during the resistive evolution, so the current sheets have a finite width rather than 
being true discontinuities.) 

The top row of figure [ID] shows / and f 2 for the initial condition of the new 
simulation. For computational reasons, we use a rectangular domain — 6 < x < 6, 
— 6 < y < 6 and —24 < z < 24. The initial background field ( 1211 has been scaled in 
r and z so that the cylindrical flux tube boundary is located at R = 5.22 (shown by 
the circle) and the overall twist remains it. Because the superimposed twist regions are 
localised well within the flux tube, the field line mapping at r = R follows that of the 
background field (a uniform twist of tt with a small perturbation). Hence the overall 
initial field for r < R retains the property of ( T2T]) that T ; * t = and T? nt = 3. The 
bottom row of figure [10] shows that this topological property persists throughout the 
dynamical evolution, verifying that both T^ nt and T^ t are conserved, acting as dynamical 
constraints. The simulation is imperfect in that the cylindrical flux tube boundary does 
not remain perfectly ideal, due to the finite numerical dissipation and to the fact that 
it does not coincide with the computational boundary. Thus the period 2 points on the 
flux tube boundary move slightly between the initial and final states. Nevertheless, it 
remains sufficiently close to ideal to illustrate the topological constraints. 
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(a) Initial state 




Figure 10. Colour maps for the first two iterations of / in (a) the initial state and (b) 
the final state of the new simulation described in section 14.21 This is a magnetic field 
with two independent dynamical constraints. Outside the flux tube R = 5.22 (shown 
by the black circle), white points indicate field lines which leave the side boundary of 
the computational box (where they are "line-tied" ) . 
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Appendix A. Proof of Lemma 13.31 

This relies on the boundary condition (JT&]) . and on the property that / maps the 
boundary 8Dq to itself. Without loss of generality we work in Cartesian coordinates 
(x, y) where the x-axis corresponds to dD and the upper half of the xy-plane to the 
interior of D . The fixed point is at x — y — 0. Let Vf be the local linearisation of Vf 
about the fixed point, so that V/ = Ax, where the matrix A is given by 

A= ( a b A — Df(0, 0) - Id. (A.l) 
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Figure Al. Possible local behaviour of a generic boundary fixed point, where Ai < 0, 
showing (a) a "semi-saddle" point and (b) a "semi-node" . The latter is possible only 
if | Df ^ 1 at the fixed point. 



We start by noting that, since / maps the boundary y — to itself, (1,0) T must 
be an eigenvector of A. This implies that a = Ai and c = 0, where Ai is the eigenvalue 
corresponding to the direction along the boundary, and that both eigenvalues are real. 
It follows from the relation Tr A = Ai + A2 that d = A2. 

Next, we derive a constraint on A2 from the boundary condition ( Ff6|) . This condition 
may be expressed as f y (x, y) > for y > 0. Since f y (x, 0) = 0, for small y we may write 

df y 

f y ~ 

dy 

Hence we must have d/V^/l(o,o) > 0. Forming the product of both sides of (lA.ip with 
x = (x,y) T then gives the following equality for the ^-component: 

-l)y = X 2 y, (A.3) 



y. (A.2) 

(0,0) 




(0,0) 

from which we conclude that 

A 2 > -1. (A.4) 

Finally, we use the assumption that \Df(0, 0)| = 1. Writing fj,i, /i 2 for the 
eigenvalues of Df (0,0), this gives /xi/i 2 = 1. The eigenvalues of Df (0,0) are related to 
those of A by Ai = \i\ — 1, A 2 = yU 2 — 1, so that 

A 2 = ^ (A.5) 

If Ai > 0, then clearly A 2 < 0. If Ai < 0, then we can use (1A.4I) to see that 1 + Ai > 0, and 
so A 2 > 0. Thus Ai, A 2 are real and of opposite sign. This rules out the possibility of a 
"semi-node" on the boundary (figure [AT]) . There remain two types of generic boundary 
fixed point: both "semi-saddles" but differing in the sign of Ai. The two types are shown 
in figure El along with the form of the extended vector field in each case. It is evident 
from the figure that one type has index while the other has index —1. □ 
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